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Force  missions.  In  this  study  modern  control  theory  was  applied  to  determine  the  time 
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that  coast,  boost,  and  sustaining  thrust  could  all  be  optimal.  For  the  particular  case 
of  sustaining  flight  it  was  shown  that  the  position  vector,  the  velocity  vector,  and  ch= 
mass  of  the  rocket  vehicle  must  satisfy  a  certain  explicitly  obtained  relationship,  and 
that  the  variable  chrust  magnitude  depends  strongly  on  the  aerodynamic  characteristics 
of  the  vehicle  and  the  optimal  trajectory  flown.  Particular  attention  was  given  to  the 
solution  for  minimum  time  lateral  turns  in  a  horizontal  plane  for  a  lifting  missile. 

A  qualitative  analysis  was  then  carried  out  for  some  typical  missions.  The  numerical 
solutions  indicate  tint  the  precentage  performance  gain  with  throttling  increases  with 
the  difficulty  of  the  mission  (especially  position  target  missions.) 
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I,  INTRODUCTION 

This  report  presents  the  results  of  a  study  to  evaluate  thrust  magni¬ 
tude  control  (TMC)  in  Air  Force  missions.  In  this  study,  modern  control 
theory  was  applied  to  determine  the  ime  history  of  the  thrust  profile  to 
achieve  optimal  maneuvers  of  aerospace  lifting  vehicles. 

This  report  consists  of  five  parts.  After  an  introduction,  we  give  a 
general  discussion  of  optimal  thrust  magnitude  control  in  Chapter  II.  A 
mathematical  model  for  a  very  general  propulsion  system,  the  system  (S) , 
was  formulated.  From  this  system  we  deduce  the  most  commonly  used 
system,  the  system  (S^ ) ,  where  the  exhaust  velocity  of  the  gas  ejected 
from  the  engine  is  constant.  Optimal  control  laws  for  the  thrust  direction, 
and  the  thrust  magnitude  are  obtained.  It  is  shown  that  coast,  boost,  and 
sustaining  thrust  all  may  be  optimal.  In  the  singular  case  of  sustaining 
flight  it  is  shown  that: 

1.  Along  the  sustaining  flight  path,  the  position  vector  T,  the  velocity 
vector  V  and  the  mass  m  of  the  rocket  vehicle  must  satisfy  a  certain  rela¬ 
tion  obtained  explicitly. 

2.  The  variable  thrust  magnitude  depends  strongly  on  the  aerody¬ 
namic  characteristics  of  the  vehicle,  and  the  optimal  trajectory  flown. 

The  results  in  this  chapter  were  obtained  for  a  Newtonian  gravitation¬ 
al  force  field.  In  the  remaining  part  of  the  report  we  consider  in  detail  the 
particular  case  of  turning  flight  in  a  horizontal  plane,  in  a  uniform  gravi¬ 
tational  field. 

In  Chapter  111  we  present  an  analytical  approach  to  the  solution  for 
minimum  time,  lateral  turns  in  a  horizontal  plane  of  a  lifting  missile.  To 
simplify  the  problem,  it  is  assumed  that  the  thrust  direction  is  always 
aligned  with  the  velocity  vector.  The  lift  and  the  bank  controls  and  the 
thrusting  program  for  the  three  types  of  optimal  arcs  involved  are  ex¬ 
pressed  in  terms  of  the  state  variables  and  a  set  of  constants  of  integration. 
Hence,  the  variational  problem  is  resolved  ultimately  to  a  multi-point 
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boundary-value  problem.  For  the  case  of  rectilinear  flight,  and  for  the 
case  of  turning  flight  with  the  final  position  being  free,  it  is  shown  that 
the  variable  tlirust  magnitude  is  given  by  (see  Table  of  Nomenclature  for 
the  notation) 


“-‘I 

Oo 

where  f  is  a  known  function  of  V/c  and  a  constant  of  integration  "a".  Fui'- 
thermore,  the  function  f  is  insensitive  to  "a".  Hence  this  formula  displays 
explicitly  the  variation  of  T  in  terms  of  the  flight  velocity,  or  equivalently, 
as  a  function  of  Mach  number.  We  notice  that  T  is  proportional  to  the  zero 
lift  drag  coefficient  ,  and  the  atmospheric  mass  density  p.  Hence,  the 
variable  thrust  magnitude  for  sustaining  flight  is  an  exponential  function  of 
the  altitude.  This  formula  gives  a  complete  criterion  lor  the  programming 
of  the  variable  thrust  profile. 

Chapter  IV  gives  the  numerical  results  of  the  same  problem.  In  this 
chapter  a  discussion  of  the  computer  program  used  to  compute  solutions  to 
fully  constrained  problems  is  presented  along  with  simulation  results  for 
representative  cases.  While  in  Chapter  III,  to  ease  the  analytical  discus¬ 
sion,  it  is  assumed  that  the  aerodynamic  characteristics  Cj^^  and  k  in  the 
parabolic  drag  polar  representation 

Cjj  =  Cjj^  (M)  +  k(M)Cj^  a.  2) 

were  independent  of  the  Mach  number,  the  numerical  program  discussed  in 
Chapter  IV  has  provision  for  inclusion  of  these  variations. 


% 


V 

% 


In  the  concluding  chapter.  Chapter  V,  we  discuss  the  five  types  of 
rocket  motors  used  in  Air  Force  missions  in  order  of  increasing  thrust 
controllability.  A  qualitative  analysis  was  carried  out  for  some  typical 
missions.  The  numerical  simulations  indicate  that  the  percentage  perform¬ 
ance  gain  with  throttling  increases  with  the  difficulty  of  the  mission  (espe¬ 
cially  position  target  missions). 
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n.  OPTIMAL  THRUST  MAGNITUDE  CONTROL 

Optimal  thrust  magnitude  control  for  flight  in  a  vacuum  has  been  dis¬ 
cussed  by  Leitmann  (Ref.  1)  and  Marec  (Ref.  2).  For  flight  inside  an  atmo" 
sphere,  partial  results  mainly  concerning  a  rocket  engine  with  constant 
exhaust  velocity  motmted  rigidly  fixed  to  the  vehicle  flying  in  a  flat  earth 
model  have  been  obtained  by  several  authors  (e.g. ,  Bryson  and  Lele  (Ref. 
3)). 

In  this  chapter  we  extend  the  results  of  Marec  concerning  a  general 
propulsion  system  to  the  case  of  flight  in  a  general  force  field  and  in  a  re¬ 
sisting  medium. 

II.  1  The  Propulsion  Systems 

Let  T  be  the  thrust  developed  by  the  engine.  It  is  assumed  that  the 
direction  of  the  thrust  can  be  taken  arbitrarily.  It  will  be  shown  that  in  this 
case  the  optimal  thrust  direction  is  parallel  to  the  vector  p^,  associated  to 
the  velocity  vector.  The  two  .remaining  control  parameters  are  the  thrust 

dull 

magnitude  T  and  the  mass  flow  p  ^  . 

Let  us  define  a  general  propulsion  system  (S)  in  the  fr,P)  space  (Fig, 
n,  1).  Let  c  be  the  exhaust  velocity  of  the  gas  ejected  from  the  engine. 

The  thrust  magnitude  is  then 

T  =  pc  01, 1) 

The  power  cf  the  propulsion  system  is 


P=4(-m)c' 


at.  2) 


The  power  being  limited  by 


Pi  P 


max 


ai.3) 


the  control  domain  in  the  (T,p)  space  is  bounded  by  the  parabola  OABC 
given  by  the  equation 


T*  =  2P  p 
max'^ 


ai.4) 


In  the  general  system,  system  (S) ,  we  assume  that  the  exhaust  velocity  is 


3 


f^Hiucii»i  lii . ■III,..- .  . . 


bounded  by 


c  .  i  c  ?i  c 


ai.5) 


This  further  restricts  the  control  domain  in  the  (T,p)  space.  In  general 
the  mass  flow  is  bounded  by 


P  =  P 
*^max 

and  the  resulting  control  space  is  shown  in  Fig,  II.  1. 

We  notice  that  the  thrust  magnitude  is  bounded  (Fig.H.  1) 

T  ^  T  =  min(T„,T  ) 
max  B  c 


(II.  6) 


(11.7) 


From  the  general  propulsion  system  (S)  we  have  the  following  special  cases: 

If  the  exhaust  velocity  is  constant,  c  .  =  c  =  c,  we  have  the  sub- 

mm  max 

system  (Sj )  (J'ig,  II.  2),  In  this  system,  the  only  control  parameter  left  be¬ 
sides  the  thrust  direction  is  either  the  mass  flow  p,  or  the  thrust  magnitude 
T.  The  thrust  is  bounded  by 

OSTiT  ai.8) 


For  an  ideal  electric  propulsion  system,  the  subsystem  (Sj)  (Fig. II.  3), 


we  have 


0,  c  =  ,  P  -  « 

max  '^max 


(n.9) 


The  control  space  is  bounded  by  the  limitation  of  the  power  only. 


In  the  following,  we  shall  examine  the  thrust  magnitude  control  of 
these  models  with  a  detailed  analysis  of  the  system  (Sj ) ,  namely  a  rocket 
engine  with  constant  exhaust  velocity. 

U.  2  Optimal  Thrust  Magnitude  Control 


The  motion  of  a  vehicle,  considered  as  a  mass  point,  flying  in  a  gen¬ 
eral  gravitational  field  and  subject  to  aerodynamic  and  thrusting  forces,  is 
governed  by  the  equations 

dir 

dt 


I  i 
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^  =~(T  +A)  +  g(ir,t)  (n.io) 

dm  _  T 
dt  “  '  c 

Using  the  maximum  principle,  for  a  minimizing  problem,  we  form  the 
Hamiltonian 

H  =  p_..V-l-p^.(i(T+A)+il-p^'2  (II.U) 

where  p^,  p^  and  p^  are  the  adjoints  to  r  ,  V  and  m. 

For  optimal  thrust  control,  we  maximize  H  with  respect  to  T.  Con- 

sider  the  vector  p  ,  called  the  primer  vector.  For  the  direction  of  T, 

maximizing  H  is  equivalent  to  selecting  T  such  that  the  dot  product  p^  •  T, 

is  maximized.  Hence  we  have  the  generalization  of  Lawden's  law  for  flight 

in  a  vacuum,  namely  that  the  direction  of  T  should  be  orientated  along  the 

vector  p  and 
V 

max(p^  •  T)  =  p^T  (U.  12) 

where  p^  is  the  magnitude  of  T.  Using  this  condition  and  upon  replacing  c 
by  T/p  we  rewrite  the  Hamiltonian  (II.  11) 

Pv 

H=-PmP  +  ^T  +  G  (U.13) 

where  G  is  independent  of  the  control  elements  p  and  T. 

We  shall  consider  successively  the  maximization  of  H  in  the  control 
space  (T,p)  for  the  different  propulsion  systems  (S),  (S^)  and  (S2). 

II.  2. 1  Optinaal  Thrust  Magnitude  Control  Using  the  System  (S) 

At  each  instant  t,  the  equation  H  =  constant  is  the  equation  of  a  straif^t 

mp 

line  in  the  plane  (T,P)  (Fig.  11.4).  The  slope  of  this  strai^t  line  is  -r — . 

Fv 

The  operating  point  in  the  (T,p)  space  depends  on  the  value  of  this  slope. 
When  p  =  0 


Dsssasww! 


m 


P=o  p^ 


(H  -  G) 


(II.  14) 


Since  m/p^  is  positive^  to  xnaximize  H  we  select  the  operating  point  p  in 
the  bounded  domain  (T,p)  representing  the  system  (S)  such  that  T  is 


maximum  possible.  Hence  we  have  the  following  cases. 


p=o 


mp 


1.  If 


m 


flight). 


2.  If 


mp 


m 


>  c^^^  we  should  select  the  point  O.  T  =  0  (coasting 


^max'  operating  point  can  oe  anywhere  in  the  seg¬ 


ment  OA.  We  have  T  =  variable,  with  c  =  c  .  We  have  the 

max 


case  of  sustaining  fli^t.  (Singular  case) 


c  mp 

,  „  max  ^  *^m 

3.  If  — = —  1  — 


<  noticing  that 


roex 


is  the  slope  of 


the  tangent  to  the  parabola  AB  at  the  point  A,  we  deduce  that  the 
maximum  of  H  is  obtained  by  using  the  point  A.  The  thrust  mag¬ 
nitude  is  constant  (Boosting  flight  at  constant  thrust) 


2P 


T  =  T.  = 
A 


max 


max 


The  corresponding  mass  flow  is 


2P 


P  = 


max 


max 


(11.15) 


(U.16) 


C„  C  mp 

.  .  B  c,  ^  m  ^ 

4.  If  max(-T~,  -r-)  ^ 2 
Z  2  D 


max 


,  the  operating  point  is  on  the 


arc  AB  (or  AC)  of  the  parabola  (II.  4).  The  thrust  roagnitude  is 
variable  and  is  given  by 


T  = 


P  P 
max*^v 

mp 


(H.  17) 


Hence  it  is  determined  by  the  ratio  p^/mp^^.  The  corresponding 
xnass  flow  is 


P 


^  *^m 


max 


(n.18) 


I 


5.  If 


.  <  max{— 5“  the  operating  point  is  at  the  point  B  (or 

P  £t  Q 

C).  The  thrust  magnitude  is  at  tlie  znaximum  constant  thrust 
(Boosting  flight  at  maximum  thrust). 


Hence  for  a  system  (S)  we  have  five  types  of  optimal  arcs. 


1 .  T  =  0  (Coasting  Arc) 

2.  T  =  variable  at  low  level  Sustaining  Arc  at  low  thrust) 

3.  T  =  constant  at  low  level  (Boosting  Arc  at  constant  thrust) 

4.  T  =  variable  at  high  level  Sustaining  Arc  at  high  thrust) 

5.  T  =  (Boosting  Arc  at  maximum  thrust) 

Besides  the  singular  case  (case  2)  which  cannot  be  treated  by  first  order 
theory,  the  thrust  magnitude  control  is  governed  by  the  quantity  mp^/p^. 


II. 2. 2  Optimal  Thrust  Magnitude  Control 


The  optimal  thrust  magnitude  control  for  the  system  Sz )  can  be 
easily  deduced  from  the  results  concerning  the  more  general  system  S)* 
The  operating  point  is  always  on  the  parabola  given  by  Eq.  (II.  4),  Hence 
the  thrust  magnitude  can  be  taken  as  the  control  available,  the  engine  being 
operated  at  maximum  power  for  any  mass  flow.  The  optimal  thrust  mag¬ 


nitude  is  given  by 


V 

T  =  — ^  P 
mp  max 


(11.19) 


Hence  it  is  a  function  of  the  quantity  p^/mp^. 


U.  2. 3  Optimal  Thrust  Magnitude  Control  Usint 


Finally,  let  us  consider  the  propulsion  system  represented  by  the 
model  (S| ) .  The  exhaust  velocity  being  constant,  the  thrust  magnitude  can 
be  used  as  the  unique  control.  Using  the  same  geometric  approach  as  for 
the  system  (S)  we  can  easily  show  the  three  types  of  optimal  arc 

mPjjj 

— =>  c,  we  take  T  =  0  (coasting  arc) 

^v 


1.  If 


2.  If 


c,  we  take  T  =  variable  (sustaining  arc).  This  case  is 


the  singular  case  and  the  thrust  magnitude  cannot  be  decided  by  the 

first  order  theory, 

mp 

3.  If -  <  c,  we  take  T  =  T  (boosting  arc). 

Py  max 

Analytically «  using  the  optimal  law  (II.  12)  for  the  thrust  direction,  we  re¬ 
write  (II.  11)  for  this  case  where  c  =  constant. 


H=p  ‘V+p  *g  +  — to  ‘A)  + — (p 
*^r  *^v  ®  m^v  m  V 


ai.20) 


We  define  the  switching  function 

mp, 

K  =  p - : 

Then,  to  maximize  H  with  respect  to  T: 


ai.21) 


if  K  <  0,  we  select  T  =  0 

if  K  =  0,  we  select  T  =  variable 

if  K>  0,  we  select  T  =  T 

max 

We  notice  that,  if  the  gravitational  field  is  time  invariant,  the  Hamiltonian 
H  is  a  constant  of  the  motion. 

To  derive  the  expression  for  the  variable  thrust  magnitude  control, 
we  notice  that  in  this  case  K  ==  0,  and  we  have  the  relation 


Py"  — 


(U.22) 


At  this  point,  it  is  necessary  to  make  assumptions  concerning  the  atmo< 
spheric  force  A  and  the  acceleration  of  the  gravitational  field  g. 


We  assume  that 


A  ='ipSV*a 


where  p  is  the  atmospheric  mass  density 

o  -  n  ^-Mr-R) 
P  -  Poe 


(n.23) 


ai.24) 
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where  R  is  a  reference  length,  usually  taken  as  the  radius  of  the  planet  and 
X  is  the  hei^t  scale  constant.  We  shall  consider  a  Newtonian  gravitational 


field 


2  =  -Ji.  £ 

g  r2  r 


(n.25) 


With  these  assumptions  we  liave  the  equations  for  the  adjoints 

*^^r  _  jt. r**  _  -«►,  ,  ^7*-  £ 

dt  ”  ?  m  ^  r 

dt  "  "  Pr  '  mV^  ^ 


ai.26) 


1  r*-  ^1.1  rp 

" .  ■■■  =  —~2  (p  •  A)  +  “~5  p  T 

Now,  along  a  sustaining  arc,  the  equation  (11.22)  is  identically  satisfied  for 
a  finite  time  interval.  Hence  we  can  take  its  derivative  to  have 

=  2!^  +  !“  42=  J_s  .I,  +  Xto 


dt  c  dt  c  dt  me  '^v 


me'^v  c 


Using  (II.  22)  we  have 


1 

“dT  ^  mT’ 


ai.27) 


On  the  other  hand 


Py  dt  "  Py  *  dt  ”  Pr  ’  Py 


2(Pj.  •  A)(p^  •  V) 


(n.28) 


By  eliminating  dp^/dt  between  the  last  two  equations,  we  have  another  re¬ 
lation  for  sustaining  flight 


(p  .  A)  2c(p  •  V) 

^  -r  '*^v  r  ^  '^^V  T 

Pr  ‘  Py  me  ^v  V*  ^ 


ai.29) 


Using  the  Hamiltonian  (11.20)  along  a  sustaining  arc  we  can  give  this  rela¬ 
tion  another  form  independent  of  the  aerodynamic  force  ^ 

^  2c{p  .V) 

-cSj.  •  p^  =  [H  -  Pj^ .  V+  ^5^  •  r)3[p^+ - ^ - ] 


ai.30) 


If  the  adjoints  and  p^  can  be  found,  this  equation  gives  the  relationship 
among  the  state  variables  Ir,  V  and  m  along  a  sustaining  arc.  Special  cases 
of  this  relation,  mainly  for  flat  earth  model  and  horizontal  flight  have  been 
obtained  by  various  authors  (e.g. ,  Ref.  3). 

We  also  notice  that,  if  the  aerodynamic  force  is  vanishingly  small, 
then  from  (11.29) 


p  •  p  =0 


ai.31} 


We  have  the  classical  result  for  orbital  transfer  in  a  -^cuum,  namely  along 
a  sustaining  arc  where  the  thrust  magnitude  is  variable,  the  two  adjoint 
vectors  p^  and  p^  are  orthogonal. 

The  thrust  magnitude  along  a  sustaining  arc  is  obtained  by  taking  the 
derivative  of  the  singular  relation  (11.30).  After  some  algebraic  manipula¬ 


tion  we  have 


(11.32) 


We  see  that  if  (p^»A)  ^  0,  the  variable  thrust  magnitude  can  be  obtained 

from  this  relation.  Some  special  cases  of  this  equation,  giving  T  explicitly 
in  terms  of  the  dimensionless  velocity  V/c  will  be  presented  in  the  next 
chapter. 
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Fig.  II.  1  Propulsion  System  (S) 


Pig.  II.  2  Propulsion  System  (Si) 
(c  =  constant) 


Fig.  II.  3  Propulsion  System  (SJ 


Pig.n.4  Optimal  Thrust  Magnitude  Control  for  System  (S) 


m.  ANALYTICAL  SOLUTIONS 


In  this  chapter  we  present  an  aimlytical  appr9ach  to  the  solution  for 
minimum  time,  lateral  turns  in  a  horizontal  plane  of  a  lifting  missile.  The 
lift  and  the  bank  controls  and  the  thrust.ng  program  for  the  three  types  of 
optimal  arcs  involved  will  be  expressed  in  terms  of  the  state  variables  and 
a  set  of  constants  of  integration.  Bence,  the  variational  problem  is  resolved 
ultimately  to  a  multi-point  boundary-value  problem. 

Simple  cases  where  certain  end- conditions  are  relaxed  can  be  solved 
analytically.  Furthermore,  in  these  cases,  the  variable  thrust  pi'ofiles  are 
obtained  expli'sitly^  thus  allowing  constructive  suggestions  for  the  design  of 
TMC  rockets. 


III.  1 .  Equations  of  Flight  in  a  Borizontal  Plane 

The  equations  of  motion  for  a  coordinated  turn  in  a  horizontal  plane, 
with  the  thrust  always  aligned  with  the  velocity,  are  (Fig.  Ill  1) 


dt' 


V  cosp 


dt> 


V  sinp 


iY  =  -  PI 

dt’  m 


=  P 

dt’  m 


dm  _  T 
dt’  "  ■  C 


aii.i) 


The  notation  is  given  in  the  nomenclature  section.  Note  that  we  use  t’  to 
designate  the  real  time.  For  the  flight  to  stay  in  the  horizontal  plane,  we 
have  the  constraining  relation 

L  coso-  =  mg  (III.  2) 

The  drag  polar  used  is  a  parabolic  drag  polar 

S  =  an.  3) 
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R 


t 


IM 


k-: 
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In  order  to  obtain  an  analytical  solution,  we  shall  assume  that  the  zero-lift 

drag  coefiticient  and  the  induced  drag  factor  k  are  independent  of  the 
Do 

Mach  number  end  the  Reynolds  number.  This  is  especially  true  for  flight 
in  the  hypervelocity  regime.  Aerodynamic  configurations  and  flight  regimes 
where  these  coei^icients  depend  on  the  Mach  number  will  be  considered  in 
the  next  chapter  where  a  complete  numerical  solution  for  a  specified  type  of 
missile  will  be  presented.  We  shall  use  the  usual  assumption  for  the  lift 
and  drag  forces  of  the  form 


an.  4) 


L  =  -ipSC,V2 

J.< 

D=^pSCj.V= 

To  obtain  general  optimal  laws  for  a  general  type  of  missile  we  intro¬ 
duce  the  nondimensional  quantities 


W  *= 


“S- 


T  = 


•fpS(?C 


Do 


011.5) 


V 

*^=7 ' 


^•s/kC 


Do 


1 

2E* 


If  the  angle-of-attack,  or  equivalently  the  lift  coefficient  C^,  is  not  con¬ 
strained,  a  natural  choice  for  the  aerodynamic  control  would  be  the  bank 
angle  v.  Then,  the  lift  and  the  drag  coefficients  will  be  obtained  from  (Ul,2) 
and  011.3).  In  the  practical  case  where  the  lift  coefficient  is  bounded,  the 
corresponding  bound  on  the  bank  angle  is  a  function  of  the  state  variables 
m  and  V  through  the  relation  Off*  2).  Hence,  the  lift  coefficient  will  be  a 
better  choice  as  control  parameter  in  this  case. 

We  define  a  lift  control  parameter  k  such  that  is  the  ratio  of  the 
induced  drag  to  the  zero-lift  drag. 


kC* 

X*  = 


on.  6) 


Then  we  l*ve 


*1 

■a 


•4 

I 


I 
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(HI.  7) 


_  w 


Xu* 


cos  or  =  ^  ,  Sintr  =  -  ■  -r.:y”  * 


=  u* 

w 


We  notice  that  when  X  =  1 ,  the  flight  is  effectuated  at  maximum  lift-to-drag 
ratio  E*. 

The  independent  controle  are  t  and  X.  They  are  bounded  in  the  space 


OStSt^ 

osxsx^ 

The  state  equations  become,  with  the  nondimensional  variables, 

dx 


an.  8) 


dt 


=  u  cos  P 


*  =usinp 

^  a+xV) 

^  _  n/X^  U*  -  W^ 
dt  uw 

dw 
dT  " 

^  =1 
dt 


an.  9) 


To  write  the  variational  equations  for  optimal  trajectories,  we  introduce 
the  adjoint  components  p^,  (i  =  0, . . .  ,5),  to  form  the  Hamiltonian 


H  ^  p.  ^  u(p.oos  p  ^p,8inp)  -  -nR  ^ 


•  /w* 


on.  10) 

The  Pj,  (i  =  1, . . .  ,5),  are  respectively  paired  with  the  first  five  equations 
OU.  9)  while  po  is  associated  to  the  last  equation.  These  adjoint  variables 
satisfy  the  equations 


dt 


=  0 


16 


dPi 

dt 


d.'p? 

dt 


dt 

dt 


0 

0 


.  .  2ti{1  +  X.2)u  (X.^u^+w^) 

-(P.cosp+p,smp)*— 


u(pjsinp-  P2COS  p) 

2}I^  _ 


an.  11) 


w‘ 


Pi  + 


w 


P4 


The  solution  to  the  problem  is  obtained  by  integrating  the  system  of  state 
equations  an.  9),  and  adjoint  equations  an.  11)  with  tlie  appropriate  end- 
conditions,  while  selecting  t  and  X.,  subject  to  the  constraints  an.  8),  in 
such  a  way  that,  at  each  instant,  the  Hamiltonian  H  given  by  aH.lO)  is  an 
absolute  maximum. 


111.2.  Optimal  Controls 

The  problem,  as  formulated,  has  a  number  of  integrals  of  the  motion. 
First,  H  does  not  contain  t  explicitly.  Since  the  final  time  is  not  specified 

K  s  constant  =  0  an.  12) 


we  see  that  Po  =  constant,  and  for  minimum  time  problem 

Po  =  -  1 

Also ,  from  the  adjoint  equations 

p,  =  aj  =  constant 


Pj  =  aj  =  constant 

Using  these  relations  we  write  the  equation  for  04 


dt 


dy  dx  d  ,  . 


dt  dt 


Upon  integrating 


P4  =a,y-a2X  +  a4 


011.13) 


an.  14) 


an.  15) 


where  a4  is  a  constant.  We  notice  that  the  constants  obtained  thus  far  are 


valid  over  the  entire  optimal  trajectory 

By  inspection  of  the  Hamiltonian^  we  define  the  switching  function 


K  =  Pj  -  wps 


(in.  16) 


Then,  for  the  thrust  magnitude  control,  to  maximize  H  with  respect  to  t: 


an.  17) 


K  >  0,  we  select  t  =  Tj^  (boosting  phase) 

K  <  0,  we  select  t  =  0  (coasting  phase) 

K  =  0,  for  a  finite  time  interval, 

we  select  t  =  variable  (sustaining  phase) 


For  the  lift  control,  considered  as  function  of  X,  H  reaches  a  maximum 


either  at 


X  =  X, 


an.  18) 


or  at  an  interior  point  given  by  ^  =  0.  Explicitly,  we  have 


(in.  19) 


Hence,  if  the  maximum  lift  coefficient  is  used,  we  have  for  the  bank,  the 
maximum  possible  bank  angle  for  horizontal  flight 


cos<r  = 


X  u'' 
M 


an.  20) 


In  the  case  where  variable  lift  coefficient  is  used,  the  bank  angle  is  given  by 


tan  <r  ■-  -  St—. 
2tip3  w 


011.21) 


We  shall  refer  to  the  interior  lift  and  bank  as  the  normal  lift  and  normal 
bank  controls.  The  other  optimal  possibilities  will  be  referred  to  as  the 
maximum  lift  and  the  maximum  bank  controls.  We  shall  assume  that  there 
are  no  wavering  maneuvers,  the  heading  being  continuously  increasing  along 
a  turning  flight,  and  hence  o-  >  0.  This  in  turn  implies  tliat  pj  and  P4  are 
of  the  same  sign. 


We  shall  analyze  separately  each  of  the  three  types  of  optiro  .  arcs 
III-2-1  Coasting  Arc 


A  =  u*  -  wr' 


{ni.22) 


A  =  aicosp  +  a2Sinp 

Along  a  coasting  arc,  t  =  0,  and  we  write  the  Hamiltonian 


uA  =  1  + 


011.23) 


The  lift  control  is  either 


X  =  X-.  or  A  =  " 

M  2iiR, 


011.24) 


We  first  consider  the  case  where  ^  The  equations  of  motion,  with  u 

as  the  new  independent  variable  are 

w  =  constant 


w  cos  0 

tiO+Xm)  u 

w  sinp 

qO+xj^j)  u 


011.25) 


^  -  M 
du  ■  "  •n0+x|^)u^ 


TiO+xyu^ 


The  boundary  conditions  at  the  ends  of  the  coasting  arc  are 
t  =  to,  u  =  uo,  X  =  xo,  y  =  yo,  0  =  Po*  w  =  Wo 
t  =  t,,  u  =  Ui,  X  =  x„  y  =  yj,  P  =  Pi,  w  =  wp 
By  integrating  the  last  equation  OU.  25)  we  have  for  the  velocity 


on.  26) 


’Xl  +  X*  )  '  u  uo 


011.27) 


The  change  in  the  velocity  from  the  initial  point  is  obtained  by  replacing 
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w  =  constant 


sTC^^k  77T^’k>^ 


1 


t 

e 


I 

I 


I 

! 


f 


dg  _  uA _ 

du  ti[A^  +  (u^+  w*)] 


it 

du 


wu^ _ 

tj[A^  +  (u*+  w^)] 


Since  P4  is  given  in  terms  of  x  and  y  by  (111.15),  Eq.  (III.  32)  can  be  solved 
for  A  in  terms  of  the  state  variables  x,y,u  and  p.  Subsequently  the  first 
three  equations  (III.  33)  can  be  integrated  numerically.  By  satisfying  the 
terminal  conditions  at  the  ends  of  the  coasting  arc  ohe  can  deduce  the  values 
of  the  constants  a^  ,82  and  a4  and  the  three  new  constants  of  integrations  b^ , 
b2  and  b4  arising  from  the  integration.  The  minimum  time  of  flight  is  ob¬ 
tained  by  integrating  the  last  equation  (111.33). 

Hence,  for  a  pure  coasting  flight,  when  the  final  states  x,y,u  and  p 
are  all  constrained,  the  optimal  flight  is  at  variable  angle  of  attack.  On  the 
other  hand,  when  these  values  are  free,  the  flight  is  effected  at  maximum 
angle  of  attack. 

We  notice  that  if  the  optimal  trajectory  also  includes  other  fypes  of 
arcs,  then  the  constants  of  integration  aj  ,82  and  84  remain  the  same 
throughout  the  entire  trajectory,  while  the  constants  of  integration  b^ ,  b2 
and  b4  are  functions  of  the  values  of  the  state  variables  at  the  ends  of  the 
coasting  arc.  Their  determination  will  specify  the  switching  point  joining 
a  coasting  arc  and  a  thrusting  arc. 

III-2-2  Sustaining  Arc 

Along  a  sustaining  arc,  K  =  0,  and  we  have 

R,  -  wps  =0  (III.  34) 

Consequently,  relation  (III.  23)  is  also  valid  along  this  type  of  arc. 


A 


i 


24 

1 


We  first  consider  the  case  where  this  arc  is  flown  with  maximum  lift. 
Consider,  in  general 


21 


W 


W 


wA 


Along  a  coasting  arc,  t  =  0,  and  a  sustaining  arc,  K  =  0,  we  have 
d(wp5)  , 

dt  w  "  wA 

On  the  other  hand,  from  the  adjoint  equations  (III.  11) 


an.  35) 


dpa  A  .  2ii(l  +  X^) 

=  -  A  +  — - - -  UPj  -  - - y-r—  P4 

dt  w  ”  wu^A 


an.  36) 

By  taking  the  derivative  of  an. 34)  and  using  (III.  35)  and  an.  36)  we  have 

A  =  p,  -  an.  37) 


W 


By  eliminating  P)  between  an. 23)  and  an.  37),  and  putting  \  =  't.^,  we  have 
a  relation  among  the  state  variables  along  a  sustaining  arc  flown  at  maxi¬ 
mum  lift 

uwA[(u  +  2)  -  u(u+l)A]  =  (ajy- a2X+a4)[A*  -  w*(u  +2)]  (in.38) 

Upon  taking  the  derivative  of  this  equation,  using  the  state  equations  for 
simplification  we  can  express  the  variable  thrust  magnitude  control  in 
terms  of  the  state  variables.  Explicit  expressions  for  the  thrust  profile  in 
some  special  cases  of  interest  wiU  be  presented  in  the  next  section. 


For  the  case  where  the  sustaining  arc  is  flown  with  variable  lift  co¬ 
efficient,  we  use  the  second  relation  ail.24)  to  rewrite  ail.23)  and  ail.37) 
we  have 


uA 


=  1  - 


[A^  -  (u*  ^w^)] 
2uwA 


P4 


an.  39) 


and 


A  = 


-A^ 

2uwA  u 


OH.  40) 


The  relation  among  the  state  variables  along  a  sustaining  arc  using  variable 
lift  is  obtained  by  eliminating  the  lift  control  A  between  these  two  equations. 
We  have  for  the  vari.?.ble  lift  control 


Hence,  if  impulsive  thrust  is  permitted,  we  have,  along  an  impulsive -thrust 
arc 


we“  =  C 


an.  48) 


111.3.  Solutions  for  some  Special  Cases 


We  have  shown  in  the  preceding  section  that  for  the  general  case,  the 
optimal  controls  can  be  expressed  explicitly  in  terms  of  the  state  variables 
and  a  certain  number  of  constants  of  integration.  In  other  words,  the  set  of 
adjoint  equations,  except  for  the  case  of  maximum  thrust  with  variable  lift 
coefficient,  can  be  integrated  completely.  Thus  the  variational  problem  is 
reduced  to  a  multi- point  boundary- value  problem.  Because  of  the  non¬ 
linearity  of  the  state  equations ,  the  actual  optimal  trajectory  can  be  obtained 
only  by  using  numerical  integration. 

In  this  section  we  shall  relax  some  final  conditions  on  the  state  variables 
in  order  to  get  deeper  into  the  problem  analyticaUy.  In  particular  we  are 
able  to  obtain  the  thrust  profile  along  a  sustaming  arc  in  these  cases.  This 
result,  coupled  with  the  nmnerical  mvestigation  in  the  next  chapter  will 
allow  us  to  formulate  constructive  suggestions  for  the  design  of  TMC  rock¬ 
ets. 


III-3-1  Rectilinear  Flight 

This  problem  has  been  analyzed  numerically  in  Ref. 4,  for  a  very 
special  type  of  rocket.  In  this  section  we  shall  give  the  general  analytical 
solution  to  the  problem  for  a  general  type  of  lifting  rocket. 


For  this  case  y  =  0,  p  =  0,<r=0,  and  the  lift  coefficient  is  always 

X  =  ^  i  an.49) 


The  lift  coefficient  is  no  longer  an  independent  control,  and  we  should 


an.  56) 
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which,  for  the  cases  of  coasting  flight,  and  sustaining  flight,  is  reduced  to 


Next,  we  consider 


dt  w  wu^  P’ 


(III.  57) 


an.  58) 


For  coasting  and  for  sustaining  flight,  we  have 

d(wpg)  _  11(0*  -  w^^ 


an.  59) 


We  consider  succesively  the  tbUowing  three  types  of  optimal  arc. 
Coasting  Arc 


We  have,  t  =  0,  and  hence 


w  =  constant 


an.  60) 


The  adjoint  variable  pj  is  obtained  in  terms  of  u  from  an.  57),  For  Ps  ,  we 


rewrite  an,  58) 


dps  _  11  (u*  -w^)_  _  _  (u*  -w^)(aiU‘l) 

dt  ^u^  P*  ~  w(tt*”+w*) 


an.  61) 


The  equation  for  u,  with  t  =  0  is 

du  _  n(u*  •*'W^) 


an.  62) 


Hence 


^  ^u-(atu-.l)(u*.“wfl 

du  11  (u*  +  w  ) 


an.  63) 


Since  w  is  constant,  ps  is  obtained  by  quadrature  and  the  system  of  adjoint 
equations  is  completely  integrable. 

The  velocity  is  obtained  by  integrating  (Ul.  62)  from  the  initial  point 
of  the  coasting  arc.  We  have 

n  .  1  ,  u^-  ''^2w  u  +  w  4.  1  .  .  'J2w' u  j**®  „„  ^ 

-ir(t-to)  =  r^S^iog  +  r^S^Arc  tan  (  an. 64) 


WT  u^+<v/Zwu  +  w  2^/2' 
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For  the  range,  using  u  as  vtiriable,  vte  have 


dx  _ 
du 


+  w*) 


Upon  integrating 


Sustaining  Arc 


— •  (x  -Xfl)  =  log  --5~r — 2 
w  “  u  -v  w 


(III.  65) 


Along  a  sustaining  arc,  K  =  0  and  we  have 

Rj  -  wps  =  0 


{UI.66) 


By  taking  the  derivative  of  this  equation,  and  using  (111.59)  and  the  second 
of  tlie  Eqs  (III.  52)  we  have 


*  wu’  " 


(m.67) 


By  eliminating  pj  between  (III.  57)  and  (III.  67)  we  have  the  relation  between 
the  nondimens ional  weight  and  the  lusndimensional  velocity  along  a  sustain¬ 


ing  arc. 


u[u*(u  4-  1)  -  w^(u  +  3)]  _  1 
(u  +  2)  (u*  -  ^  ai 


(in.  68) 


In  Fig.  III.  2,  we  have  plotted  the  ti’ajectoriea  for  sustaining  flight,  in  the 
(w,u)  space,  usings  =  l/a^  as  parameter.  Explicitly,  we  have  the  equation 
of  the  curves  with  a  <  0 


'/u*  +  (3-a)u-2a 


(III.  69) 


The  family  of  curves  is  bounded  by  two  limiting  curves.  When  -a  -♦  w  we 
have  the  parabola 

w  =  u^  (III.  69a) 


When  a  -*■  0,  we  have 


i  /u  •*'  1 

^  -  “Vm 


an.  69b) 


This  limiting  curve  corresponds  ix>  the  singular  curve  obtained  in  Ref.  5,  for 
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the  problem  of  maximizing  the  range. 

To  have  the  thrust  profile  along  a  sustaining  arc^  we  take  the  time 
derivative  of  (III.  69)  using  the  appropriate  relations  for  du/dt  and  dw/dt. 
Upon  simplification  we  have 


^  ,  2uVt(2:_a,)u,:2a]  ^ 

[u*  +  p  -a)u-  2a] 

f2u*  +  (9  -  4a)u^  +2(a^  -  8a4-3)u^ 
[u®  +2 (3  - a)u*  +  (a^  -  12a  +  12)u^  +6 (a*  - 


an.  70) 


-2a(7-4a>i  +  8a^] 

4a  +  l)u^  +  2a  (6a  -  7)u  +  8a^] 


In  Fig.  III.  3,  we  have  plotted  the  thrust  magnitude  of  sustaining  flight  as  a  func¬ 
tion  of  the  velocity,  using  "a"  as  the  parameter.  The  family  of  curves  is 
bounded  by  two  limiting  curves.  When  a  **“00  we  have 


T 


u  +  2 


(III.  70a) 


When  a  0 ,  we  have 


(u  +  3)(u*  +  6u*  +12U  +  6) 


an.  70b) 


From  Fig.ni.2,  it  is  seen  that,  along  a  sustaining  arc,  as  w  decreases, 
the  velocity  decreases.  From  Fig. III. 3,  we  see  that  the  variable  thrust  is 
also  decreasing  along  a  sustaining  arc.  The  level  of  the  thrust  is  high  for 
large  velocity,  and  remains  at  lower  level  for  smaU  velocity.  In  particu¬ 
lar,  for  the  range  of  velocity  studied  in  Ref. 4,  between  M  =  0.3  and  M  = 
2.5,  one  easily  verifies  in  our  plot  for  the  thrust  profile  that  t  remains 
small  and  is  slowly  decreasing.  It  has  been  found  numericaUy  in  Ref.  4 
that  the  thrust  magnitude  along  a  sustaining  arc  is  nearly  constant  and  re¬ 
mains  small. 

We  notice  that  the  exact  solution  for  our  thrust  profile  applies  to  a 

general  type  of  missile.  By  the  definition  au.  5)  of  t,  we  see  that  the  real 

thrust  is  proportional  to  (pSc*C„  ).  Hence  our  formula  (III.  70)  gives  ex- 

Uo 

plicitly  the  variable  thrust  profile  in  terms  of  the  altitude  (throu^  p ) ,  the 
flight  velocity  V,  and  the  missile  characteristics  S  and  C-.  . 


Using  the  expression  (III.  70)  for  the  variable  thrust,  we  have  the  equa¬ 
tion  for  the  distance  travelled  along  a  sustaining  arc 


.  r  2(u  +  2)ND+2u(u2  +  2a)  , 

TMlu— p-I,u-btT 


(in.  71) 


where 


N  =  u^  +  (1-  a)u-  2a 
D  =  u^  +  (3  -a)u-  2a 


011.72) 


It  is  known  that  the  integral  in  Eq.  0^*7 1)  can  be  expressed  in  terms  of 
elliptic  functions.  The  time  is  finally  obtained  from 


-4q(t-to)=  r  du 

J  un/Nd  [u2+ (2-a)u-2a] 


011.73) 


.  uo 


Hence,  the  equations  along  a  sustaining  arc  are  completely  integrable. 


Along  a  boosting  arc,  where  t  =  Tj^,  the  variation  in  the  nondimen- 


sional  weight  is  simply 


w-wo  =  -‘TT^O-to) 


011.74) 


In  the  (w,u)  space,  the  trajectory  for  boosting  arc  is  obtained  by  into- 
grafcuig  the  equation 


au.75) 


Although  an  analytical  solution  to  this  equation  is  not  known,  it  is  seen 
that  for  Tj^  -♦  » ,  the  limiting  curve  is 


u 

we  =  C 


where  C  is  a  constant  of  integration. 


Optimal  Trajectory 


The  problem  of  minimum  time  between  two  given  end  points  in  the 
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(w,u,x)  space  is  hence  obtained  by  joining  the  different  types  of  optimal  arcs. 

Let  B  •  Boosting  arc,  C  =  Coasting  arc  and  S  *  Sustaining  arc.  Follow¬ 
ing  the  same  type  of  discussion  as  in  Ref.  6,  we  see  that  there  are  six  pos¬ 
sible  types  of  optimal  trajectory,  namely  (Fig.  III.  4) 


BC 

CB 

BSC 

and 

CSB 

BSB 

CSC 

Although  the  actual  computation  involves  numerical  solution,  it  is  possible 
to  have  the  following  general  properties  which  facilitate  the  computation 
(Fig.  III.  5). 

Proposition  1 

Let  X  =  (w,  u,x). 

If  Xj,  and  Xj  are  both  specified,  tiien  in  general  the  optimal  trajectory  in¬ 
volves  all  three  types  of  arc. 

Proposition  2 

If  Wo  >  Uo 

the  initial  arc  is  a  boosting  arc 

the  initial  arc  is  a  coasting  arc. 

Proposition  3 

If  w^>  uj 

the  final  arc  is  a  coasting  arc. 

If  Wj<  u* 

the  final  arc  is  a  boosting  arc, 

III-3-2  Turning  Flight  with  Free  Lnd- point 

As  in  the  case  for  rectilinear  flight,  the  turning  flight  problem  re¬ 
duces  to  a  problem  containmg  one  single  parameter  if  we  do  not  specify  the 


end  position  x^  and  y^. 


In  this  case,  from  the  general  results  obtained  in  Section  I1I-2  we 


have 


and  hence 


aj  =  0,  a2  =  0 


p4  =  a4  =  constant 


(III.  77) 


(m.78) 


The  controls  involved  are  the  lift  (and  hence  the  bank  angle)  and  the  thrust 
magnitude.  For  the  optimal  lift,  we  have  either 


X  = 


an.  7  9) 


011.80) 

2nP3 

We  shall  consider  the  case  of  maximum  lift  and  next,  the  case  of  variable 


Maximum  Lift  Program 


The  lift  is  given  by  OU.  79).  and  hence,  we  have  for  the  bank  angle 


cos  <r  - 


X_.u' 

M 


on.  81) 


The  optimal  arc  is  of  three  types . 


Coasting  Arc 


The  general  solution  in  section  III-2  is  valid  and  we  have  the  following 
results  starting  from  the  initial  point  Wq,  Uq,  Po>  sq. 


w  =  Wo 

1  1 

- ^(t-to  ):---  — 

2na+X.\,)(P-  Po)  =  >>.log  ^ - - 

M  M  +'s/X^^u*  -  w*  “  Uo 


M  /  ^  ,  Uo 

- (s  -  So)  =  log 

Wo  u 


on.  82) 
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T  3u^(u  +  4) 

(T+X^)  "  (2ui  +  9u  +  i2) 


an.  85b) 


Once  the  thrust  is  known,  the  variation  in  the  heading,  and  the  time  of  flight 
along  a  sustaining  arc  are  obtained  by  quadratures.  We  observe  that,  along 
a  sustaining  arc,  both  the  velocity  and  the  variable  thrust  level  are  decreas¬ 
ing. 


Boosting  Arc 

When  X  =  Xj^,  the  variation  of  the  mass  in  terms  of  the  velocity  is 
given  by  Eq.  an. 47)  written  in  terms  of  the  initial  conditions  at  the  start 
of  the  boosting  arc  as 


The  variation  of  w  as  function  of  time  is  simply 

w  -  Wo  =  ‘  an.  87) 


From  these  we  have  u  as  function  of  t.  For  the  change  in  the  heading,  we 
write 


^  = 
du 


u[tj^  -  a+%)u2] 


ail  .88) 


w^  is  obtained  from  an.  86)  and  subsequently  we  have  p  in  terms  of  u  by 
quadrature. 

Normal  Lift  Program 


The  variable  lift  is  given  by  (111.80).  The  optimal  arc  using  a  variable 
lift  coefficient  is  of  three  types. 

Coasting  Arc 


Since  A  =  0,  P4  =  a*  we  rewrite  the  Eqs.  an. 22)  -  ail. 24) 
A  =  VX^u*  -  w* 


I  .  0 

w  uw 


an.  89) 
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By  eliminating  Ps  among  the  equations ,  we  have  the  equation  for  A. 

a4A*  -  2uwA  -  a4  (u^  +  w^)  =  0  (HI.  90) 


Solving  for  this  equation,  we  have  explicitly  for  the  lift  coefficient 

=  1  +  ^  N/u^w^  +  ai^  (u^  +  w*)  (III.  9 1) 

For  the  bank  angle  we  have  (a4  f  0) 


tan  a 


u‘ 


+  a^(l  + 


an.  92) 


Since  w  is  a  constant  along  a  coasting  arc,  we  can  plot,  for  different  values 
of  w,  the  lift  coefficient  X ,  and  the  bank  angle  cr,  in  terms  of  the  flight  ve¬ 
locity  u,  using  a4  as  parameter.  For  a  given  missile,  w  can  be  taken  as  the 
initial  weight  or  the  final  weight.  To  have  the  plots  valid  for  a  general  type 
of  missile,  it  is  convenient  to  include  the  constant  a4  in  the  variables.  Let 

U  .  ^  ,  W  .  ^  au.93) 

Then  we  have  the  equation  for  X  . 

X^  =  1  +  a  +  u*)  ±  +w‘  +  u‘w‘  ail.  94) 

u«  '  /  —  u-* 


with  the  condition 


i  <V 

u* 


M 


an.  95) 


We  can  have  a  general  plot  for  the  lift  coefficient,  in  terms  of  the  velocity 
using  w  as  parameter. 

Similarly,  using  u  and  w  we  have  for  the  bank  angle  along  a  coasting 
arc  using  variable  lift  coefficient 


tan  a 


u  ± 


on.  96) 


Since  w  is  constant,  we  see  from  this  formula  that  when  u  -*•  0  <r  -*•  45*  and 
when  il  -*  eo,  v  90“ .  Since  TI  is  decreasing  along  a  coasting  arc,  the  bank 
angle  for  normal  lift  is  also  decreasing.  By  £q.  (III.  8 1)  we  see  that  the 
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bank  angle  is  also  decreasing  along  a  coasting  arc  when  maximum  lift  co¬ 
efficient  is  used.  The  analytical  results  obtained  in  this  chapter  are  in 
agreement  with  the  results  from  an  independent  numerical  analysis  carried 
out  in  Chapter  IV. 

Sustaining  Arc 

For  a  finite  time  interval  of  flight  along  a  sustaining  arc,  we  have 


p,  -  wps  =0 

The  lift  coefficient  is  given  by  (111.42)  with  A  =  0.  We  have 


(m.97) 


an.  98) 


On  the  other  hand,  from  Eq,  HI. 44,  by  putting  A  =  0,  P4  =  a4  ,  we  have  the 
relation  between  the  weight  and  the  velocity  along  a  sustaining  arc  using 


variable  lift 


u^w^(u+  2)(u*  -  w^)  _2 


on.  99) 


Explicitly,  using  a  =  a^  as  parameter,  we  have 


w  =  u 


2  /[u*(u+2)^♦^2a(u^M)]•^u^^/at  +  2)tu^(u•^2)  +  4a 
V  2[u3(a  +  2)  +  a(u+l)2] 


an.  100) 


Fig.  HI.  8  plots  the  family  of  sustaining  arc  in  the  (w,  u)  space.  The  curves 
are  bounded  by  two  limiting  curves. 

For  a  -•  0,  we  have  the  parabola 


w  =  u* 


For  a  00,  we  have  the  curve 


an.  lOOa) 


on.  100b) 


Using  an.  100)  we  can  express  the  variable  lift  control  in  terms  of  u,  and 
the  constant  a  as 


=  1  + 


u^(u+2)-t-2a(u+l)  -  u^y(u^^■:^)Tu^fu^•2)^■4a 
[u*(u  +  2)+a(u+l)*3 


an.  101) 
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For  u  -*•  0,  X  ^/3'.  For  u  -*  oo,  X  -♦  1.  Hence,  we  see  that,  when  variable 
lift  is  used  along  a  sustaining  arc 

l<X<^/3  (in.  102) 

The  family  of  lift  curves  is  bounded  by  two  limiting  curves.  When  a  -♦  0, 

Eq.  (HI.  101)  becomes 

x  =  i  an.ioia) 

When  a  -*  00 ,  we  have  the  limiting  curve 


X 


[u  -H  3 
V  U  +  1 


an.  101b) 


For  the  bank  angle,  we  have  in  general 


Using  ail. 98)  we  have 


tan^o-  =  -^(X^u*  -  w^) 


W.  =  filiSi  (^  - 1) 

u  W^ 


(m.io3) 

an.  104) 


Another  expression  for  tan<r  can  be  obtained  by  noticing  that  Eq.  ail.  90)  is 
also  valid  for  sustaining  flight,  and  that  A  =  w  tana.  Hence  we  can  write 
the  equation  for  tana  with  a4  ^  0 

tan^a  -  ~  tana  -  a  +  -^  )  =  0  ail.105) 

Combining  the  last  two  equations  we  have 

tana  =  -  (u  +  1)]  011.106) 

For  a  to  have  positive  value,  we  notice  that  a4  <  0,  when  w  >  u^/'s/u+T . 

By  eliminating  w^  between  the  two  equations  011.104)  and  011.106)  we  have 
the  equation  for  tana. 

tan^o-  -  iilHl^tana-  (u  +  2)  =  0  011.107) 

a4 

Since  a4  <  0,  we  take  the  positive  root 


tana  =  -^^(0+2) -n/(u+2)[u‘(u+2) +  4aJ]  3 


011.108) 


Fig.  in.  9  plots  the  variation  of  the  bank  angle  along  a  sustaining  arc 
in  terms  of  the  nondimensional  velocity,  using  a4  as  parameter. 
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When  u-»  0,  tan  o-  "»  •s/T and  <r  -*  54"44* . 

When  u‘^00  (r'*90“or  0®  depending  on  the  sign  of  a4 . 

The  family  of  curves  is  bounded  by  the  limiting  curve  obtained  when  -a4  *♦  «>. 
We  have 

tan  <r  =  •sJvT+'Z  (III.  108a) 

The  variable  thrust  profile  is  obtained  by  taking  the  derivative  of  the  Eq. 

(Ill,  99),  using  the  appropriate  relations  for  du/dt  and  dw/dt.  After  some 
algebraic  manipulation,  we  have 

^  ^  2uf(u  +  l)D-N  ^  uH(2u  +  5)N-D]  +  a(3u  +  4)f(u  +  l)N-D] _ 

D  u3[(2u2  +  6u+5)N- (u  +  l)2D]+a(2u2+5u+4)[(u  +  l)N- D] 

(m.l09) 


where 


N.  =  [  u^  (u  +  2)  +  2a  (u+ 1)]+  u^  n/(u  +  2)[u*(u  +  2)  +  4a  ] 
D  =  2[u’(a+2)  +  a(u+ 1)?] 


(III.  110) 


Fig.  III.  10  plots  the  variation  of  the  variable  dimensionless  thrust  t  in 
terms  of  the  dimensionless  velocity  u,  using  a  as  parameter.  The  family  of 
curves  is  bounded  by  two  limiting  curves. 

When  a  •*  0  .  , 

(ni.l09a) 


T  = 


u  +  2 


When  a  00 


2u^(u-f2}(3u-<-4) 
^  (u  +  l)(2u2 +5u  V4) 


(III.  109b) 


It  is  seen  that  t  is  increasing  with  u.  Hence,  since  along  a  sustaining  arc, 
the  velocity  is  decreasing,  the  variable  thrust  is  decreasing  along  that  arc. 
Also  we  notice  that  along  the  sustaining  arc,  for  small  velocity  the  variable 
thrust  magnitude  remains  at  relatively  low  level. 

Boosting  Arc 

For  variable  lift  control,  no  analytical  solution  can  be  found  for 
the  boosting  arc,  but  it  is  seen  that  in  the  (w,  u)  space  the  trajectory  tends  to 
the  exponential  curve  (III.  48)  when  Tj^  -♦  oo. 


Optimal  Trajectory 

As  for  the  case  of  rectilinear  flight,  the  optimal  trajectory  is  of 
the  six  .types 

BC  CB 

BSC  and  CSB 
BSB  CSC 

The  optimal  trajectory  is  obtained  by  joining  different  types  of  subarcs  such 
that  the  end  conditions  are  satisfied. 

To  fecilitate  the  computational  program,  we  may  use  the  following 
criteria  (Pig.lU.ll). 

Proposition  1 

When  maximum  lift  is  used 


If  w.  >  (LjjUj 
The  initial  arc  is  a  boosting  arc. 


If  Wo  < 


7uo  +  3 


the  initial  arc  is  a  coasting  arc. 
When  normal  lift  is  used 
If  Wo  >  u^ 

the  initial  arc  is  a  boosting  arc. 

the  initial  arc  is  a  coasting  arc. 
Proposition  2 

When  maximum  lift  is  used 

»  W{  > 

the  final  arc  is  a  coasting  arc. 

If  w.<  -JSLL 


U  0 


C  B 

Hjr.m.4 


o  CSB 


CSC 


Optimal  Trajectories  for  Rectilinear  Flight 
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Fig.  Tii.  6  Optimal  Turn  Trajectories  for  Sustaining  Flight  at 
Maximum  Lift 


Fig.  III.  7  Variable  Thrust  Profile  Along  Sustaining  Arc  for  Turning 
at  Maximum  Lift 


Fig.  III.  8  Optimal  Turn  Trajectories  for  Sustaining  Flight  with 
Variable  Lift 


Maximum  lift 


Normal  Lift 

Fig.  m.  J 1  Criteria  for  the  Selection  of  the  Initial  and  Fim.l  Arcs 
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IV.  NUMERICAL  RESULTS 


In  this  cliapter  a  discussion  of  the  computer  program  used  to  compute 
solutions  to  fully  constrained  problems  will  be  presented  along  with  simula¬ 
tion  results  for  representative  cases. 

IV.  1.  Basic  Optimization  Problem 

The  equations  of  motion  for  flight  in  a  horizontal  plane  are  given  by 
Eqs.  (III.l)  and  (III. 2).  To  facilitate  the  discussion  of  this  section,  the 
state  variables  (i.  e. ,  x,  y,  v,  p,  m)  will  be  denoted  by  Xj , . . . ,  x  ,  where 
n  =  5  for  this  problem,  and  the  controls  T.<r,Cj^  will  be  denoted  by  Uj,U2,U3,  re¬ 
spectively.  The  computer  program  is  developed  in  the  (x,  u>-system. 


The  problems  considered  in  this  study  are  Mayer  type  optimal  con¬ 
trol  problems,  and  the  computer  program  uas  developed  for  the  following 
Mayer  problem. 


MINIMIZE: 

J  <|>(tj,Xj) 

av.i) 

SUBJECT  TO: 

/s/ 

X  =  f(t,x;Ui  ,U2,U3) 

av.2) 

x(to)  =  Xo 

av.3) 

0  i  u,  S  T 

*  max 

av.4) 

U3  sc 

^max 

av.5) 

G(X,U2,U3)=0 

av.6) 

Equation  (IV.  1)  is  the  performance  index  which  contains  the  quantity  to  be 
optimized  (e.  g. ,  t^  for  minimum  time  problems)  and  penalty  terms  Tor  the 
terminal  boundary  conditions.  Equations  (iV.  2),  (IV. 4)  and  (IV,  5),  and 
(IV. 6)  correspond  to  Equations  (III.l),  (111.3),  and  (III. 2),  respectively. 
The  control  equality  constraint  (IV.  6)  can  be  used  to  eliminate  either  Uj  or 
U3  from  the  problem.  Since  the  main  purpose  of  the  program  is  to  verify 
and  extend  the  results  of  Chapters  II  and  III,  Cj^  is  eliminated,  i.e. , 

G(X,U2,U3)=0  =?>  U3=g(X,U2),  (IV. 7) 

or,  in  physical  variables- 
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r---w.  < - 


L  cos  or  =  mg  => 


'L  p  S  cos  (T 


av.8) 


We  now  have  the  following  problem;  Minimize  Eq.  (IV.  1)  subject  to 
x(to)  =  xo  and: 

X  =  fft.x^Uj.Uz)  ^"fLt.x^Ui.Uz.g^x.Uz)]  (IV. 9) 


0  2  Uj  s 

g(x.  Uz)  i  C, 


(IV.  10) 
(IV.  11) 


IV.  2 .  Solution  by  Gradient- Type  Methods 

The  first  step  in  the  development  of  a  function  space  gradient- type 
method  for  optimal  control  problems  is  the  definition  of  the  following  aug¬ 
mented  functional: 


1 

J  =  <!>(t^.x^)  +  J  \'^(t)[f(t,x,Ui,U2)  -x]dt 


av.12) 


where  X(t)  is  a  vector  of  influence  functions  which  in  the  limit  (of  the  numer 
ical  algorithm)  approaches  the  p(t)“vector  of  Chapter  III. 

It  is  convenient  to  define  a  Hamiltonian  function 

H(t,x,\,u)  =\^f(t,x.u).  av.l3) 

Let  u^**^(t)  a  (Ui°\t),U2*^(t))  be  an  estimated  control  vector  which  satisfies 

(0) 

Inequalities  (IV.  10)  and  (IV.  1,1),  and  let  t^  '  be  an  estixxiate  of  the  final  time . 
One  can  then  integrate  Eqs.  (IV.  9)  forward  to  define  a  corresponding  tra¬ 
jectory  x^°^(i),  defined  on  [to,ti°^],  and  a  value  for  the  perforaaance  index. 

(0) 

Let  \  '(t)  be  an  arbitrary  continuous  vector  which  we  shall  character¬ 
ize  later.  To  develop  the  gradient- class  of  numerical  algorithms,  expand 
Sq.  (IV.12)  to  first-order  about  the  pair  (x^**\t)  ,u^^\t))  with  the  following  defi¬ 


nitions 


u4t)  =  u^°H)  +  6u{t) 
x^^H)  =  +  6x(t) 


aV,14) 


aV.15) 


I  Xj  I 

.  jtf  -  Ak]it  . 

to 


av.ii 


Upon  integration  by  parts  of  the  third  term  in  the  integrand  and  rearrange 
ment  of  terms,  the  following  equation  expresses  the  difference  in  cost  be¬ 
tween  the  base  trajectory  and  the  first  iterate  (to  first-order): 


AJ[6u]  =  +H^®^)^^^jdtj  +  (<j,j^|^-X^®^iJ^^dx. 


»  f 


J^^x  H^®^'^6u3dt  . 


(IV.  IS 


where 


6x(t’®^  =  dXj  -  x^^(t^^)dtj 


(IV.  2C 


*(o)^ 


has  been  used  to  eliminate  6x(t^')  from  Eq,  (IV.  19). 

,(o). 


We  now  characterize  X  '(t)  so  that  a  stable  iterative  algorithm  (i.  e. , 
AJ  i  0)  is  defined. 


SPECIFY: 


=  -H Jt.x%),u^®V),X^®^t)]  . 


av.z; 

av.22 


Then, 


AJ  =  (/^  +  H^®Ldt,  +  Ch!®^  6udt, 


(IV.  2: 
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I  T 

=  J{[  -  to)]dtj  +  6u}dt.  av.24) 

to  ^ 

Note  that  if  the  following  choices  are  made: 

dtf  =  t^'^  -  tj")  =  +  H^“h/(tJ^-  to)]  (IV.  25) 

6u(t)  =  u^^^(t)  -  u^°^(t)  =  -€H^®^(t)  (subject  to  constraints) 

av.  26) 

with  €  >  0  (small),  then  AJ  ^  0.  These  choices  represent  the  gradient 
method  choices. 

In  the  computer  program  developed  here,  the  conjugate  gradient 
method  [Refs  7,8]  was  employed.  It  can  be  shown  that  this  method  also 
guarantees  AJ  ^  0.  We  shall  now  present  the  algorithm  in  two  parts;  first, 
the  algorithm  with  no  control  inequality  constraints  will  be  listed,  and  sec¬ 
ond,  the  modifications  necessary  to  handle  the  control  inequality  constraints 
will  be  listed 

UNCONSTRAINED  CONJUGATE  GRADIENT  ALGORITHM 

1.  )  Guess  u^°\t) ,  t^'  . 

2. )  Compute:  x^\t),  X^^(t),  H^\t) 

t^  T 
^.ffla) 

p<I)(t)  =  H^^t)  +  - P^'^^t)  aV.27) 

j.fjj(I-i)  jj(I-1)jH 

/  u  u 
to 

p^“^(t)  s  H^^t) 
ta) 

3. )  Check dt  s  € .  If  yes,  stop.  If  no,  goto  4). 

to 

4.  )  Perform  one- dimensional  (1-D)  search  to  determine  the  value aj  which 

minimizes 

J[u^^  -  ajp^^]. 


av.28) 
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5.)  Define  set  1  =  1+1.  Return  to  2.) 

Equation  (IV.  27)  defines  the  conjugate  gradient  search  direction.  Dis¬ 
cussions  of  this  method  are  given  in  Refs  [7,8]. 

MODIFICATION  FOR.  0  ^  Ui  i  T 
_ *  max 

Let  be  the  set  of  points  at  which  Ui^^(t)  is  on  the  boundary,  i.  e. , 


W,**’  2  {tlu,%)  =  0  or  . 

Replace  the  Ui  -part  of  the  inner  products  in  Eq.  (IV.  27)  by 


f  H<'>H<'>dt. 

J  “l  “l 


av.29) 


av.30) 


i.e. ,  in  the  evaluation  of  the  inner  product  associated  with  H  ,  (t)  do  not 

“1 

include  boundary  subarcs. 

The  set  is  formed  during  the  1-D  search.  The  implementation 

is  as  follows.  Let  a  be  a  candidate  search  parameter; 


if  u ^\t)  -  a p^^  >  T  .set  (t)  =  T 


max 


if  u^^t)  -  ap^^<  0  ,  set  u^'^^^(t)  =  0. 


av.3i) 


After  the  function  evaluations  for  the  1-D  search  are  computed,  a  value 
is  obtained  by  a  cubic  interpolation  process.  The  set  is  then  defined 

as  the  set  of  all  t  which  cause  the  inequalities  in  Eq.  (IV.  3 1)  with  Oj  re¬ 
placing  a. 

MODIFICATION  FOR  gfx.Uj)  s  C- 
_ ^max 

As  noted  in  Chapter  III  the  constraint  on  bank  angle  is  a  function  of 
state  variables.  Thus,  this  control  constraint  is  treated  in  a  manner 
slightly  different  than  the  thrust,  Uj .  Let  W2  ^  be  the  set  of  points  at  which 
g(x^^(t),  U2^^(t))  =  .  In  the  computer  program  the  control  vector  is 

stored  at  each  integration  step,  and  linear  interpolation  is  used  to  define 
the  control  between  the  fixed  step-lengths.  As  with  is  deter¬ 

mined  in  the  1-D  search.  However,  in  the  forward  integration  for  the 
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function  evaluation  for  a  specified  search  parameter,  V ,  the  control  con¬ 
straint  relation  is  changing  because  x  appears  in  g(x,U2)  =  ^Lmax' 
when  the  integrator  is  at  t^,  the  constraint  at  t^^  +  At  is  predicted  by 

g[x(t  )+i(t.)  At.  U2  (t  +  At)]  ^  C,  .  aV.32) 

*  1  A 


If  the  integration  stepsize  is  sufficiently  small,  this  approximation  is  suf¬ 
ficient  to  approximate  the  control  constraint  at  t  +  At.  Except  for  this 
modification,  the  set  is  determined  in  the  same  way  that  is 

determined. 


IV.  3.  Deck  Description 


The  program  is  developed  based  on  the  conjugate- gradient  algorithm. 
A  Runge-Kutta  fourth  order  integration  scheme  is  used  to  perform  both  for¬ 
ward  and  backward  integration.  A  switch  is  imposed  such  that  it  can  be  set 
to  the  gradient  algorithm. 

The  program  is  designed  to  minimize  a  weighted  performance  index 
which  includes  the  following  effects. 


1.  Total  flight  time  t^  . 

2.  Terminal  states  x^.y^.v^.p  ^  and  m^. 

The  performance  index  is 

J  =  ctj+Pi(x-Xj)’  +  P2(y-yj)^+P3(v- Vj)2 +P4(p  -  Pj)2 +p5(ni-mj)^ 


SEKALF 
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MAIN:  Reads  in  all  necessary  data,  initial  and  terminal  values,  controls 
the  application  of  the  conjugate  gradient  algorithm,  calls  the  forward  and 
backward  integration  routines,  directs  the  one  dimensional  search,  updates 
the  control  vector  and  terminal  time  and  prints  out  a  message  concerning  the 
results  of  the  iteration  and  prints  out  the  control  profiles  obtained  by  that 
iteration. 

A.  J  =  An  integer  which  indicates  where  the  data  is  located. 

B.  Namelist  input  data 

G  =  gravitational  constant 

RHO  =  atmospheric  density 

CDO  =  zero-lift  drag  coefficient  (Cq^  ) 

CE  =  exhaust  velocity 

DK  =  induced-drag  coefficient  (k) 

AREA  =  aerodynamic  reference  area 
BELTS  =  integration  stepsize 
CLA  =  slope  of  the  lift  curve 

ITMAX  =  limit  on  number  of  conjugate  gradient  iterations 
ITMX  =  limit  on  steps  in  1-D  search 
KOUNTM  =  limit  on  iterations  for  wei^t  cutoff 
IKEY  =  number  of  iterations  to  reset  the  search  direction 
CSTR  =  estimate  of  final  cost,  initial  guess  of  parameter  a  for  one¬ 
dimensional  search 
PFUN  (5)  =  penalty  coefficient  vector 
CCOST(l)  =  coefficient  in  cost  functional 
B  =  not  used 

I  =  Index  for  output  device 

JERK  =  switch  for  scheme  optloo  I  » 

DTFM  =  maximum  allowable  final  time  change 
XDTFM  =  fraction  of  DTFM  used  to  start  1-D  search 
XO  =  initial  position  in  x- direction 
YO  =  initial  position  in  y- direction 
UO  =  initial  velocity 
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BO  =  initial  heading 

WO  =  initial  mass 

TO  =  initial  time 

XF  "  final  position  in  x-direction 

YF  =  final  position  in  y- direction 

UF  =  final  velocity 

BF  =  final  heading 

WF  =  final  mass 

lOUT  =  print  frequency  for  forward  integration 
lOUT  2  =  print  frequency  for  backward  integration 

IPRNT  1=1  for  initial  run,  =  2  for  continuing  run  after  normal  termina¬ 
tion  of  conditional  cutoff 

IPRNT  2  =  not  used 

CBND(2,2)  =  control  bound 

ERRMX  =  error  tolerance  for  integration  routine 
ERRMN  =  not  used 

TCUT  =  upper  time  limit  on  trajectory 

EPST  =  cutoff  tolerance  for  norm  of  control  change 

EPSTF  =  not  used 

EPSA  =  cutoff  tolerance  for  integration  weight  cutoff 
EPSIT  =  cutoff  tolerance  on  gradient  norm 
ERR  =  cutoff  tolerance  for  small  cost  change 
Note:  for  any  final  state  unspecified  an  arbitrary  number  other  than  zero 
may  be  assigned. 

C.  Control  Vector  Date 

IJKU  =  total  number  of  control  points 
TF  =  initial  guessed  t^ 

U  (3 ,  UKU)  =  control  vector  and  time  point 

D.  Bounded  Thrust  Value 

NN  =  number  of  points  with  bounded  thrust 
TAUMAX(N)  =  maximum  thrust  limit  at  Nth  point 
TAUMIN  (N)  =  minimum  thrust  limit  at  Nth  point 
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SEiCALF:  (One- dimensional  search  subroutine):  Determines  the  parameter 
for  the  new  control  value  in  the  conjugate  gradient  algorithm.  Fits  a  cubic 
in  a  to  known  values  of  J(o'),  9J/9a,  to  obtain  min  J(«)  and  then  a*  for  J  min. 

FWDINT:  Subroutine  performs  the  forward  integration  of  the  state  variables 
and  calls  the  subroutines  to  evaluate  the  cost  functional  and  final  multiplier 
values. 

BAKINT:-  Subroutine  performs  the  backward  integration  of  the  state  variables 
and  multiplier  equations ,  calls  on  GRADFN  to  calculate  the  gradient  and 
store  the  value  at  each  integration  step,  determines  the  new  search  direction. 

DERIV 1:  Subroutine  calculates  the  time  derivatives  of  the  state  variables 

DERIV2:  Subroutine  calculates  the  time  derivatives  of  the  multipliers 

AERODN;  Calculates  aerodynamic  parameters 

XLAMFN:  Calculates  final  multiplier  values 

COSTFN:  Calculates  cost  functional 

GRADFN:  Calculates  gradient  values  at  each  step. 


IV.  4.  Representative  Problems  and  Numerical  Solutions 


Several  minimum-time  missions  have  been  selected  and  solved  by  the 
numerical  programs.  The  missile  is  assumed  to  have  a  parabolic  drag 
polar 

Cj5  =  +  k(M)Cj^  (  33) 

In  terms  of  the  angle  of  attack,  the  lift  coefficient  is  given  by 


Cl  =  Cj^^(M)a 

where  Cl^  is  the  slope  of  the  lift  curve 

9C^(M,a). 

K  ~  V  A=o 


(IV.  34) 


av.35) 


The  induced  drag  coefficient  k(M)  can  be  related  to  the  aerocynam- 
ic  stability  derivatives  by  observing  that  Cjj(M,a)  can  be  expanded  in  Tay¬ 


lor's  series  as 


Cd  =  C^^(M)  + 


/®^D  \  .  1/®  ^  2 


(IV.  36) 
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For  a  parabolic  representation,  as  given  by  (IV.  33),  the  coefficient  of  the 
second  term  on  the  right  hand  side  of  £q.  (IV.  3 6)  is  negligible  and  we  have 

kmcl  =  k(M)Cl^ma^ 


Hence 


where 


Cd  W 

■  za  (M) 

Lt 

a 


av.37) 


a=0 


It  is  customary  to  write 


(M)  =  6  (M)  Cl  (M) 
aa  a 


aV.38) 


aV.39) 


where  « (M)  is  the  aerodynamic  efficiency  with  typical  values  bounded  by 

-i  ^  €  (M)  S  1  (IV.  40) 


Hence 


=  .  *  (M) — 

'  '  C.  (M) 

LIq, 


(IV.  41) 

For  a  given  Mach  number,  the  maximum  lift-to-drag  ratio  is  given  by 


E* 


1 _  =i  /  Cl,(M) 

^  ^«(M)Cdo(M) 


aV.  42) 


To  verify  the  analytical  solutions  obtained  in  Chapter  III,  in  the  first  part 
of  this  numerical  analysis.  Part  A,  we  shall  assume  that  the  coefficients 
Cdjj,  ^  constant.  We  shall  indicate  in  Part  B  the  necessary 

modification  in  the  subprograms  when  these  aerodynamics  derivatives  vary 
as  functions  of  the  Mach  number. 
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A.  Constant  Aerodynamic  Derivatives 

The  following  values  are  used  in  the  numerical  computations. 

Initial  weight.  Wo  =  861  lbs 

Final  weight,  W^  =  434  lbs 

Maximum  thrust,  Tj^  =  34,000  lbs 

Minimum  thrust,  T  ,  =  0  lbs 

mm 

Exhaust  velocity,  c  =  8,050  ft/ sec 
Reference  area,  S  =  0.66  ft^ 

Air  mass  density  at  at  40,000  ft,  p  =  0.000585  slugs/ft? 

Zero  lift  drag  coefficient,  =0.3 

Lift  curve  slope,  Cj^  =  10.3 

Induced  drag  coefficient,  k  =  0.097 

Maximum  angle  of  attack,  =  30® 


Problem  1;  Pure  Coasting  Flight 


_ Initial  Conditions _ 

Xo  =  0 

Yo  =  0 

Po  =0 

Vo  =  2136. 2  ft/sec  ~  2. 2  Mach 
Wo  =  861  lbs 


Terminal  Conditions 
Xj  =  free 
Yj  =  free 
p^  =  45® 

Vj  ^  1000  ft/sec 
W^  =  861  lbs 


This  case  was  considered  to  test  of  the  efficiency  of  the  program.  After  12 
iterations  we  obtain  =  10. 904  seconds  with  the  final  velocity  being 

V(t£)  =  1525  ft/sec.  As  predicted  by  the  analytical  solution,  the  flight  is 
effected  at  maximum  angle-of-attack,  and  the  maximum  allowable  bank 
angle  decreases  monotonically  along  the  optimal  trajectory. 

Problem  2:  Minimum  time  turning  with  high  initial  velocity 


Initial  Conditions 
Xo  =  0 
Yo  =  0 

Po  =0 


Terminal  Conditions 
Xj  =  free 


Pj  =  135* 


Vo  =  2136.2  ft/sec 
Wo  =861  lbs 


Vj  ^  1000  ft/sec 
Wj  =  434  lbs 


In  this  problem,  a  thrusting  phase  is  involved.  Since  the  constraining  final 
velocity  is  low  the  thrust  profile  is  of  the  Boost- Coast  type.  The  minimum 
time  obtained  is  t^  =  9.49  seconds.  The  trajectory  is  flown  with  maximum 
angle-of-attack. 

Keeping  the  same  final  heading,  if  we  increase  the  final  constraining 
velocity,  there  exists  a  critical  final  velocity  such  that  the  thrust  profile 
reverses  to  the  Coast-Sustain- Boost  type. 

Problem  3:  Minimum  time  turning  with  low  initial  velocity 


Initial  Conditions 
Xo  =  0 
Yo  =  0 

Po  =-0 

Vo  =  1000  ft/ sec 
Wo  =  861  lbs 


Terzxiinal  Conditions 
Xj  =  free 
Yj  =  free 
Pf  =  45“ 

V^  i  1000  ft/sec 
Wj  =  434  lbs 


After  40  iterations,  we  obtain  tf  .  =  8. 133  seconds  with  the  final  velocity 

*mm 

being  V(t£)  =  5700  ft/sec.  The  optimal  thrust  profile  is  Boost-Coast  with 
the  trajectory  flown  at  maximum  angle-of-attack. 

Problem  4:  Minimum  time  turning  to  a  specified  terminal  position 


Initial  Conditions 
Xo  =  0 
Yo  =  0 

Po  =0 

Vo  =  2136.2  ft/ sec 
Wo  =  861  lbs 


Terminal  Conditions 
Xj=  31680  ft  =  6  miles 
Yj  =  31680  f t  =  6  miles 
=  free 

Vj  i  1000  ft/sec  ■ 

Wj  =  434  lbs 


This  problem  is  designed  to  force  the  appearance  of  a  sustaining  phase 
where  variable  thrust  control  is  used.  After  23  iterations  we  obtain 
^min  “  ^  seconds.  The  trajectory  is  composed  of  an  initial  coasting 

arc  of  19  seconds,  followed  by  a  sustaining  arc  of  5  seconds  with  a  final 
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boosting  arc  of  3  seconds.  The  thrust  level  for  the  sustaining  arc  is  very 
low  and  is  nearly  constant. 

Fig.  IV.  1  shows  the  convergence  of  the  final  coordinates  and  Y^, 

while  Fig.  IV.  2  shows  the  convergence  of  the  final  weight  W^. 

Fig.  IV.  3  shows  the  optimal  thrust  profile  and  Fig.  IV.  4  presents  the 
variation  in  the  optimal  bank  angle. 


Problem  5:  Minimum  time  turning  to 
Initial  Conditions 
Xo  =  0 
¥0=0 
Po  =0 

Vo  =  2136.2  ft/sec 
Wo  =  861  lbs 


specified  terminal  position 
Terminal  Conditions 
X,  =  0 

=  105600  ft  =  20  miles 
=  free 

Vjg  1000  ft/sec 
Wj  =  434  lbs 


This  problem  is  desinged  to  have  a  longer  flying  time,  thus  making  the 

variable  thrust  control  more  prominent.  After  35  iterations  we  obtain 

tf  .  =  45.029  seconds.  The  trajectory  is  of  the  type  Coast-Sustain-Coast. 
*mm 


Fig.  IV.  5  shows  the  final  trajectory,  with  the  initial  guessed  trajectory. 
The  figure  illustrates  the  efficiency  of  the  numerical  program. 

Fig.  IV.  6  shows  the  optimal  thrust  profile  and  Fig.  IV.  7  presents  the 
variation  of  the  optimal  bank  angle. 


B.  Mach  Dependent  Aerodynamic  Parameters 


The  program  has  been  assembled  to  include  the  case  where  the  aero¬ 
dynamic  derivatives  1M)  and  k(M)  in  Eq.  (IV.  33)  are  functions  of  the  Mach 
number.  The  following  numerical  data  are  available  for  the  rocket  con¬ 
sidered. 
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Drag  Coetrxci''..xit  Cj^(o,M> 


or  {deg  ) 

0.0 

0.8 

1.0 

Mach  Number 
1.05  1.5  2.0 

2.2 

3.5 

6.0 

0 

0.25 

0.20 

0.30 

0.45 

0.40 

0.33 

0.30 

0.22 

0. 15 

5 

0.30 

0.26 

0.40 

0.50 

0.45 

0.38 

0.34 

0.27 

0. 18 

10 

0.51 

0.48 

0.70 

0.80 

0.75 

0.61 

0.60 

0.48 

0.41 

15 

1.06 

1.07 

1.30 

1  41 

1.35 

1.18 

1.14 

1.00 

0.92 

20 

2.10 

2.22 

2.60 

2.67 

2.60 

2.20 

2.10 

1.90 

1.67 

25 

3.87 

4.25 

4.60 

4.68 

4.30 

3.48 

3.30 

3.05 

2.84 

30 

- 

4.90 

5.00 

- 

5.60 

- 

4.60 

- 

- 

Lift  Coefficient  CL(a.M) 


a  (deg.) 

0.0 

Mach  Number 
1.0  1.5  2.2 

2.5 

3.5 

6.0 

0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

5 

0.67 

0.62 

0.70 

0.77 

0.80 

0.81 

0.78 

0.67 

10 

1.60 

1.54 

1.80 

1.85 

1.90 

1.91 

1.94 

1.86 

15 

2.85 

2.61 

2.90 

00 

• 

3.50 

3.58 

3.55 

3.05 

20 

4.12 

3.65 

4.00 

4.93 

5.20 

5.58 

5.27 

3.81 

25 

4.85 

4.00 

6.00 

7.50 

7.80 

7.93 

6.90 

4.70 

30 

- 

6.80 

7.00 

8.00 

8.00 

- 

- 

- 

Least  squares  is  used  to  obtain  the  values  of  C  (M)  and  k{M) 

Do 

at  each  point.  Then  polynomial  regression  curve  fitting  is  applied 
to  derive  the  expressions  for  the  functions  (M),  k(M)  aid  their  deriva¬ 
tives  9Cn  /9M,  dk/dM. 

By  the  end  of  each  boosting  phase,  the  velocity  is  usually  in  the  hyper¬ 
sonic  range.  Since  theoretically  and  k  tend  asymptotically  to  constant 
values  when  M  <»,  these  functions  are  set  constant  for  M  §  6. 0. 


seconds 


I 


seconds 


Coasting  Arc 


!0 


20 


Fig.  IV,  5  Optimai  Trajectory  with  Initial  Guessed  Trajectory 
(Problem  5) 


0,000  lbs 


V.  QUALITATIVE  ANALYSIS  AND  CONCLUSIONS 


One  of  the  objectives  of  this  study  is  to  determine  the  conditions  under 
which  a  variable  thrust  magnitude  program  is  optimal,  and  in  that  case,  its 
relationship  with  missile  parameters  and  trc  jectory.  In  the  general  case, 
such  a  relationship  was  displayed  by  Eq.  (11.32)  for  the  rocket  motor  with 
constant  exhaust  velocity.  For  horizontal,  rectilinear  flight,  the  variable 
thrust  magnitude  was  given  explicitly  in  terms  of  the  flight  velocity  in  Eq. 
(III. 70),  and  for  the  case  of  turning  flight  with  free  terminal  positions,  it  was 
obtained  throu^  the  Eq.  (III.  85)  for  flight  at  maximum  lift  and  the  Eq. 

(III.  109)  for  flight  using  variable  lift  coefficient.  In  all  these  cases,  the 
thrust  magnitude  is  always  decreasing  along  the  sustaining  flight  path  and 
furthermore,  it  remains  at  low  level  for  low  Mach  number.  These  conclu¬ 
sions  have  been  substantiated  b>  an  independent  numerical  study  carried  out 
in  Chapter  IV. 

If  an  ideal  rocket  motor  where  the  thrust  magnitude  level  is  fuUy  con¬ 
trollable  (motor  represented  by  the  system  (S^)  in  Chapter  II)  is  available, 
then  the  trajectory  can  be  flown  optimally.  In  practice,  the  rocket  motor  is 
preprogrammed  and  in  this  chapter  we  shall  compare  the  performances  of 
some  preprogrammed  motors  with  an  ideal  motor. 

We  can  classify  rocket  motors  used  in  Air  Force  missions,  in  order 
of  increasing  thrust  controllability,  as  follows  (Fig.  V.  1). 

1 .  Preprogrammed  motors.  The  thrust  profile  is  fixed  at  the  time 
of  manufacture. 

2.  Preprogrammed  pulse  motors.  The  thrust  profile  is  fixed  at  the 
time  of  manufacture;  however,  the  off  time  betv/een  pulses  may  be  con¬ 
trolled. 

3.  Stop- restart  motors.  This  system  permits  motor  shutdown  and 
restart  at  any  time  on  command.  The  thrust  level  is  either  rated  thrust  or 
zero. 

4.  Step- thrust  motors.  This  system  has  the  capability  of  providing 
thrust  at  more  than  two  discrete  levels.  The  on-off  times  and  the  order  of 
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thrust  levels  are  variable  on  command. 

5.  Continuously  variable  motor.  This  system,  idealized  by  our 
model  (Si)  offers  the  capability  to  provide  any  .thrust  time  profile  within 
specific  limits  of  thrust  and  rate  change  of  thrust. 


As  has  been  displayed  explicitly  by  the  equation  (11.32),  whenever 
variable  thrust  profile  is  optimal,  it  var.  •  '  s  function  of  the  state  variables 

r,  V  and  m  and  the  adjoint  components  p^  and  p^.  The  adjoint  p^  is  in¬ 
volved  only  in  the  switching  function  K  (Eq.  II.  21)  indicating  the  timing  for 
stop  and  restart.  It  is  kjiown  that  the  adjoint  p^  and  depend  on  the  termi¬ 
nal  conditions.  Hence  for  each  specified  mission  we  have  a  resulting  opti¬ 
mal  trajectoiry  and  a  specified  thrust  profile. 


The  numerical  examples  given  in  Chapter  IV  clearly  show  that,  unlike 
the  simple  case  of  horizontal  rectilinear  flight,  the  number  of  variables 
involved  in  the  optimization  problem  for  minimum  time  turning  flight  requires 
an  impoztant  program  for  trajectory  analysis  in  order  to  do  an  adequate 
comparative  qualitative  analysis  among  the  different  types  of  motors.  Al¬ 
though  this  task  is  laborious,  we  believe  that  with  our  analytical  results 
and  numerical  programs,  the  analysis  can  be  easily  carried  out  if  sufficient 
computational  time  is  allowed. 


In  this  concluding  chapter  we  shall  give  this  comparative  analysis  for 
two  specified  problems  given  in  Chapter  IV,  namely  problem  2  and  4. 

We  have  seen  in  problem  4  that,  by  constraining  the  final  position 
(intercept  problem),  to  meet  the  end- conditions,  the  trajectory  usually  in¬ 
cludes  a  sustaining  arc,  where  variable  thrust  is  used.  In  general  the  vari¬ 
able  thrust  level  is  low  compared  to  the  maximum  thrust.  Hence,  in  the 
rockets  used  for  the  comparison  we  shall  consider  constant  thrust  at  differ¬ 
ent  intermediate  levels  with  variable  on-off  time.  The  results  of  the  com¬ 
parative  analysis  are  given  in  the  tables  below. 


Case  1: 

Xo  =  0  .  Yo  =  0  ,  Po  =  0  ,  Vo  =  2136. 2 ft/sec.  Wo  =  861  lbs 

=  free,  =  free,  p^  =  135“,  ^  1000 ft/sec  ,  =434  lbs 

Opti.nal  Thrust  Profile  =  Boost-Coast 
Minimum  time  =9.49  seconds 

Comparison  with  Pre-prograinmed  Motors 


Motor 

Thrust  Profile  {T  in  lbs) 

Time  (sec) 

1 

Boost- Coast 

T=  34,000  0  ^ti3.14 

T  =  0  t  >  3.14 

9.49  (minimum) 

2 

Sustain- Coast 

T  =  12,500  0  5ti9 

T  =  0  t>  9 

11.2 

3 

Sustain- Coast 

T  =  6,640  0  g  t  i  6.25 

T  =  20,000  6.25  <t^9.25 

T  =  0  t>  9.25 

12.04 

4 

Sustain- Boost- Coast 

T  =  7,900  0  i  t  5  8.25 

T  =  34,000  3.25  <  t  S  9.25 

T  =  0  t>  9.25 

12. 17 

5 

.'sustain-  Boost-  Coast 

T  =  3,i30  0  i  t  i  7.  25 

T  =  34,000  7.25  <  t  i  9.25 

T  =  0  t>9.25 

13.37 

The  comparison  in  this  case  shows  the  influence  of  the  timing  and  the 
order  in  the  sequence  of  arcs. 
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Case  2: 


Xo  =  0  ,Yo  =  0  >  Po  =  0  ,  Vo  =2136.2ft/sec..  Wo  =  861  lbs 

=  31680  ft,  Yj  =  31680  ft,  p^  =  free,  ^  1000  ft/ sec  ,  =  434  lbs 

Optimal  Thrust  Profile  =  Coast-Sustain- Boost 
Minimum  Time  =  27.16  seconds 

Comparison 'With  Pre-programmed  Motors 


Motor 


Thrust  Profile  (T  in  lbs) 
Coast-Sustain- Boost 
(Fig.  IV.  3) 

Coast- Boo  St 
T  =  0  0itg24 
T  =  34,000  t>  24 
Coast-  Sustain-  Boo  st 
T=0  0Stg23 
T  =  4, 250  23  <  t  i  25 

T=  34,000  23  <tS28 


Time  (sec) 


27,16  (minimum) 


27.20 


28.157 


The  comparison  in  this  case  shows  that  the  variable  thrust  can  be 
adequately  approximated  by  either  a  null  thrust  or  a  constant  low  level  thrust. 
The  problem  of  minimum  time  turning  in  horizontal  flight  has  been  chosen 
for  our  numerical  analysis  but  the  computer  program  can  be  adapted  for 
other  types  of  performance  indices  such  as  maximum  coverage  (reachable 
sets).  Furthermore,  it  can  be  used  to  solve  optimization  problems  for 
flight  in  three-dimensional  space. 

Also  the  analytical  results  presented  in  this  report  have  been  extended 
to  the  case  of  optimal  aerodynamic  and  thrusting  maneuvers  for  three  di» 
mensional  flight  in  a  general  gravitational  force  field  (Refs.  9-11). 
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